EXACT EULER MACLAURIN FORMULAS FOR SIMPLE LATTICE POLYTOPES 
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Abstract. Euler Maclaurin formulas for a polytope express the sum of the values of a function over the 
lattice points in the polytope in terms of integrals of the function and its derivatives over faces of the 
polytope or its expansions. Exact Euler Maclaurin formulas KP2 CSl CS2 Gu BVl apply to exponential 
or polynomial functions; Euler Maclaurin formulas with remainder IKSWn IKBW?! apply to more general 
smooth functions. 

In this paper we review these results and present proofs of the exact formulas obtained by these authors, 
using elementary methods. We then use an algebraic formalism due to Cappell and Shaneson to relate the 
different formulas. 
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1. Introduction 

Let / be a polynomial in one variable. The classical Euler Maclaurin formula (see [Knl chapter XIV]) 
computes the sum of the values of / over the integer points in an interval in terms of the integral of / over 
the interval and the values of / and of its derivatives at the endpoints of the interval. The formula is almost 
three hundred years old ( \Be\ . \Eu\ . \Ma\ ) . We refer the readers to the treatments by Wirtinger |Wi| and by 
Bourbaki jBo| . 

A version of this formula that was generalized by Khovanskii and Pukhlikov to higher dimensions (see 
|KP2| ) involves variations of the interval. It reads 



(1.1) /(a) + /(a+l) + ... + /(&-l) + /(&)= Td — Td — 



d 
dhi 



d_ 

dh2 



/ll=/l2=0 



b+hi 



f{x)dx, 



-h2 



where a.b ^X, and 



(1.2) 



Td(i:>) = 



D 



1 



= 1 



-D 



12 



720 



The right hand side of (|l.l(l is well defined because /^^^^^ f[x)dx, as a function of hi and /12, is again a 
polynomial. We call this an exact formula, to distinguish it from Euler Maclaurin formulas with remainder, 
which apply to more general smooth functions. 
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In higher dimensions, one replaces the interval [a, b] by a polytope, that is, the convex hull of a finite set of 
points in a vector space, or, equivalently, a bounded finite intersection of closed half-spaces. We assume that 
our polytopes have nonempty interior. A polytope in R" is called an integral polytope, or a lattice polytope, 
if its vertices are in the lattice Z". It is called simple if exactly n edges emanate from each vertex. For 
example, a two dimensional polytope (a polygon) is always simple. A tetrahedron and a cube are simple; a 
square-based pyramid and an octahedron are not simple; see Figure^ A polytope with a non- singular fan is 




Figure 1. 

a simple polytope in which the edges emanating from each vertex lie along vectors that generate the lattice 
Z". For example, in Figure[21 the triangle on the left has a non-singular fan, and that on the right does not 
(check its top vertex). 




Figure 2. 

We refer the reader to |B1 IGrI IRWI IT) for general background on convex polytopes. 

Remark 1.1. A fan in R" is a set of convex polyhedral cones emanating from the origin, such that the 
intersection of any two cones in the set is a common face, and such that the face of any cone in the set is 
itself a cone in the set. The fan of a polytope A C M" consists of a set of cones associated to the faces of A; 
for each face we take the cone generated by the inward normals to the facets that meet at that face. 

A convex polyhedral cone is non-singular if it can be generated by a set of vectors in Z" which are part 
of a Z-basis of Z". A fan is non-singular if each cone in the fan is non-singular. The name "non-singular" 
comes from the theory of toric varieties; non-singular fans correspond to non-singular toric varieties. 

The terminology in the literature is inconsistent. A non-singular cone is also called "smooth cone" and 
"unimodular cone" . Polytopes with non-singular fans are also called "non-singular polytopes" , "smooth 
polytopes" , or "Delzant polytopes" . In our previous papers |KSW1I IKSW2| we used the terms "regular 
orthant" and "regular polytope" (not to be confused with the more common usage of this term as "platonic 
solid"). Other terms that have been suggested to us are "unimodular polytope" or "torsion- free polytope". 

Khovanskii and Pukhlikov IKPllHTF^ . following Khovanskii |KhlllKIi2] (see also Kantor and Khovanskii 
|KKl[IK^K2p . gave a formula for the sum of the values of a quasi-polynomial (polynomial times exponential) 
function on the lattice points in a lattice polytope with non-singular fan. This formula was further generalized 
to simple polytopes by Cappell and Shaneson CSl, CS2, CSS, SI, and subsequently by Guillemin Gg and 
by Brion-Vergne |BVlir5V2) . Cappell-Shaneson and Brion-Vergne jBV2| also work with polytopes that are 
not simple; in this paper we will restrict ourselves to simple polytopes. Also see the explicit formulas in jSV| 
and the survey jVelj . When applied to the constant function / = 1, these formulas compute the number of 
lattice points in a simple lattice polytope A in terms of the volumes of "expansions" of A. A sample of the 
literature on the problem of counting lattice points in convex polytopes is IMdL IVI IKKIL IMoL IPoL IDRI 
pDR, HMl IdCP] : see the survey BP and references therein. 

Remark 1.2. Khovanskii's motivation came from algebraic geometry: a lattice polytope A with non-singular 
fan determines a toric variety A/a and a holomorphic line bundle La M^. The quantization Q{AIa) 
is interpreted as the space of holomorphic sections of Xa and is computed by the Hirzebruch Riemann 
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Roch formula. The lattice points in A correspond to basis elements of Q{M/^). A simple poly tope A 
still determines a toric variety, which now may have orbifold singularities. Cappell and Shaneson derived 
their formula from their theory of characteristic classes of singular algebraic varieties. Guillemin derived his 
formula by applying the Kawasaki-Riemann-Roch formula to symplectic toric orbifolds. Brion and Vergne's 
proof uses Fourier analysis and is closer to Khovanskii and Pukhlikov's original proof. 

In this paper we present an elementary proof of the exact Euler-Maclaurin formulas that follows the lines 
of the original Khovanskii-Pukhlikov proof, through a decomposition of the polytope into an alternating sum 
of simple convex polyhedral cones. We then use an algebraic formalism due to Cappell and Shaneson to 
explain the equivalence of the different formulas. 

The proof of the exact Euler Maclaurin formula for a simple convex polyhedral cone involves the follow- 
ing ingredients: the summation of a geometric series, the change of variable formula for integration, and 
Frobenius's theorem that the average value of a non-trivial character of a finite group is zero. (See section 
01) The "polar decomposition" of the polytope into simple convex polyhedral cones was proved in papers of 
Varchenko and Lawrence We present a short direct proof of it. (See sectional) 

From a simple polytope A with d faces one gets a expanded polytope A(/i), for h — [hi, . . . , hd), by parallel 
translating the hyperplanes containing the facets, see equation H5.2|l below. The integrals of a function / 
on A(/i) and on its faces are functions of the d variables hi, . . . ,hd. The formulas of Khovanskii-Pukhlikov, 
Guillemin, and Brion- Vergne involve an application of infinite order differential operators to these functions. 
The Cappell-Shaneson formula does not involve expansions of the polytope. It is stated through a formalism 
that we call the Cappell-Shaneson algebra. Their abstract formula translates to several different concrete 
formulas; each of these involves applying differential operators to the function and integrating over faces of 
the polytope. The relations in the Cappell-Shaneson algebra allow one to pass between the different concrete 
formulas. In SectionElwe show how to incorporate expansions in h into the Cappell-Shaneson formalism. In 
section 13 we use a generalization of the "polar decomposition" , which applies to polytopes with some facets 
removed, to prove that the Cappell-Shaneson formula is equivalent to the Khovanskii-Pukhlikov formula in 
the case of polytopes with non-singular fans. 

Acknowledgement We would like to thank Itai Bar-Natan, S. Cappell, V. Guillemin, G. Kalai, G. Ku- 
perberg, and J. Shaneson for helpful discussions. We are particularly grateful to A. Khovanskii and to M. 
Vergne for explaining to us their proofs. 

2. Euler Maclaurin formulas in one dimension 

In this section we present Euler-MacLaurin formulas for a ray and for an interval, in order to illustrate 
arguments that generalize to higher dimensions. 

d 

The ODE for the exponential function. Let D — — . Since 

oh 

for any formal power series F in one variable we have 

(2.1) F{D)e^'' = F(Oe«'' 

in the ring of power series in two variables. It follows that, for any non-negative integer N, 

(2.2) F{D)^-^ = {FiOe^'^f^ 

where the superscript (N) denotes the A^th term in the Taylor expansion in ^. Under suitable convergence 
conditions, (|2.1|l is an equality of functions. See |Bo| . 

Euler Maclaurin formula for a ray. To conform with the topological literature, let us define the Todd 
function by 

(2.3) Td{S) := 

1 — e 

and the corresponding Todd operator in one variable by 

Td (D) . 
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Our general rule H2.1|l gives 

Td (D) e^^ ^ Td(Oe«'' 

in the ring of formal power series. If |^| < 2tt, so that the power series for Td(^) converges, we can regard this 
last equation as an equality of functions. Namely, the left hand side is the limit of the functions obtained by 
applying the partial sums of the infinite series Td(£') to the exponential function. If ^ 7^ 0, we can re-write 
this as 

(2.4) TdP)— = ^-^. 

If ^ > 0, the geometric series expansion 

converges, as does the integral 

h 



so (|2.4f) gives 



d 



h 



E 



(2.5) Td(^-jy e^^d. 

/!.=0 n= — 00 

This is the Euler Maclaurin formula for the ray (—00, 0], with the function f{x) — e^^ . 

Polar decomposition of an interval. In the one dimensional case, the "polar decomposition" becomes 
the relation 

(2.6) 17(0:) = 1c,(.t)- 1^9 (x) 

between the characteristic functions of the interval / = [a, 6], the ray C;, = (—00,6], and the ray = {-co, a) 
(which is obtained from the ray Cq = [a, 00) by flipping its direction and removing its vertex). 

Euler Maclaurin on finite intervals. Let / = [a, h] be a closed interval with integer endpoints. For 
h ~ {hi, ft-2), consider the expanded interval 

I{h) := [a - h2,b + hi]. 

Summation and integration of the function 

fix) = 

gives 

Mb+hi) Ma-h2) 



(2.7) I{h,i) := / e^-dx 

Ji(h) ? ? 

for all ^ such that ^ 7^ 0, and 

(2.8) S{0 E = I 



xeini 



for all ^ S C such that ^ 1. An indirect proof of (|2.7|l and (|2.8|l . which generalizes to higher dimensions, 
uses the "polar decomposition" (|2.6(l : if Rc^ > 0, then 



Ee^^dx = 7 and e^^dx = 7 = r- 



k— — oo 



Since S{^) is the difference of these two infinite sums, (|2.8() holds whenever Re^ > 0. Because the set 
G C I 7^ 1} is connected, by analytic continuation H2.8|l holds for all ^ in this set. A similar argument 
shows that (|221) holds for all f in the set G C | ^ ^ 0}. 

At this point one can proceed in several ways. 
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Formal approach. One can deduce an Eulcr Maclaurin formula for polynomial functions directly from 1)2. 7f) 
and 1)2. 8|) . This is the one dimensional case of the approach of Brion-Vergne. From (|2.7|) we get 

(2.9) eZ(/i, e) = e«(''+''i) - e«(''-''^) 
for all ^ ^ 0, and, by continuity, also for = 0. From (|2.8|) we get 

(2.10) ^SiO = Td(Oe«^ - Td{-Oe^- 

for all ^ such that ^ I and, by continuity, also for ^ = 0. Comparing the Taylor coefficients with respect 
to ^ on the left and right hand sides of 1)2.9)1 and of ()2.10)) . we get 



(2.11) 

and 
(2.12) 



N 



1(h) 



-dx 



m+h,)r+^ (e(a-M)~+i 



(iV + 1)! 



(iV + 1)! 



N 



= (e«^Td(e) - e«°Td(-e)) 



(7V+1) 



where the superscript {N + 1) denotes the summand that is homogeneous of degree + 1 in ^. Since 

d d 
Td(— — ) = 1 + a multiple of — — , 



Td 



d 
dhi 



Td 

= Td 



d 
dh2 

d 
dhi 



(Cx 



Nl 



-dx 



h=0 Jl(h) 

.,=0 (^ + 1)! 



- Td 



= (Td(Oe«^ - Td(-Oe 



d 
dh2 



^2=0 



(iV + 1)! 
by ins. 



by Em, 



Nl 



by (E321). 



This gives the Euler-Maclaurin formula 
(2.13) 



d d 



for the function f{x) = 



whenever ^ ^ 0. Because multiplication by A^! and division by the non-zero 



constant commutes with summation, with integration, and with the infinite order differential operator 

Td(g^)Td(g^), we deduce the Euler-Maclaurin formula ()2.13)) for the monomials f{x) = x^ , and hence 
for all polynomials. 

Approach through Euler Maclaurin for exponentials. In the original approach of Khovanskii-Pukhlikov, 
one deduces an Euler Maclaurin formula for polynomials, and, more generally, for (quasi)-polynomials, from 
a formula for exponentials. (A quasi-polynomial is a sum of products of exponentials by polynomials.) In 
the one dimensional case, the Euler-Maclaurin formula for exponentials asserts that 



(2.14) 

or, equivalently, that 



Td 



d 
dhi 



Td 



h=0 J 1(h) 



Td 



d 



Td 



d 
dhi 



I{h,0 = S{0. 



h=0 



This formula is true for all ^ such that |^| < 27r. For ^ in the punctured disk 

(2.15) UeC I e^O , 1^1 <27r}, 

the formula follows immediately from 1)2.6)) . 1)2. 7)1 . ()2.8I) . and from the facts that 



(2.16) 



Td 



d 
dhi 



h=0 



Mb+hi) 



^ 1 - e ^ 
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and 



(2.17) Td 



d 



(If Re^ ^ then 12.16|l is an Euler-Maclaurin formula for the ray (—00, b] or [b, 00), and similarly for l|2.17|l . 
However, and ifTTTjl hold for all ^ in the set ^J^.) 

In H2.16|l and H2.17|l . the left hand sides converge to the right hand sides uniformly in ^ on compact 
subsets of the punctured disk H2.15|l . This is because the Taylor series of Td(-) converges uniformly on 
compact subsets of the disk {|^| < 27r}, and the functions ^ and ^ are bounded on compact subsets of the 
punctured disk H2.15|l . It follows that in (|2.14|) the left hand side converges to the right hand side uniformly 
in ^ on compact subsets of (|2.15|l . But the right hand side and the partial sums of the left hand side of 
(|2.14() are analytic in ^ for all |^| < 27r. It follows from the Cauchy integral formula that the left hand side 
of H2.14|l converges to the right hand side, uniformly on compact subsets, on all of | |^| < 2tt}. 

It further follows that the infinite sum on the left hand side of (|2.14(l can be differentiated with respect 
to £, term by term. Hence, the infinite order differential operator on the left hand side of l|2.14(l commutes 
with differentiation with respect to ^. Since 

r r ^ ^ 

e^^dx— I x^e^^dx and tttt e*"^ — 'S^x'^e^^, 



we get the Euler-Maclaurin formula (|2.13|) for the function f{x) = x'^e^^ by differentiating the left and right 
hand sides of (|2.14|) k times with respect to ^. 

Approach through Euler Maclaurin for rays. Yet another approach is to deduce the Euler Maclaurin 
formula for an interval directly from the Euler Maclaurin formula for a ray. See Appendix IbI 

3. Polar decomposition of a simple polytope 

In this section we describe a decomposition of a simple polytope into simple convex polyhedral cones. 
These cones have apexes at the vertices of the polytope. Each is generated by flipping some of the edge 
vectors according to a choice of "polarization" , so that they all point roughly in the same direction, and 
removing corresponding facets. For an illustration of this decomposition in the case of a triangle, see Figure 

El 

In this section we present a short direct proof of the polar decomposition of a simple polytope, similar to 
the one that we gave in |KSW2j . In Section (see (|7.4|l 'l we give a variant of this decomposition that applies 
to a polytope with some facets removed. 



Figure 3. The polar decomposition theorem 
Let A be a polytope in an n dimensional vector space V and F a face of A. The tangent cone to A at 

is 

Cf ^ {y + r{x - y) \ r>0,yeF,xeA}. 
(Warning: other authors define the tangent cone as {r{x — y) \ r>0,y€F,xi^ ^}-) 

Example 3.1. Consider the triangle in with vertices (0,0), (2,0), and (0,1). The tangent cone at the 
hypotenuse is the (closed) half plane consisting of all points lying below the line extending the hypotenuse; 
the tangent cone at the top vertex consists of all rays subtended from this vertex and pointing in the direction 
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Figure 4. Triangle and tangent cones 



of the triangle; if F is the face consisting of the triangle itself, then the tangent cone Cp is the whole plane. 
See Figure^ 

Let ai, . . . ,(Td denote the facets (codimension one faces) of A. (Warning: Cappell and Shaneson use the 
symbols ai to denote the dual objects to the facets, namely, the one dimensional cones in the corresponding 
fan.) Assume that A is simple, so that exactly n facets intersect at each vertex. Let Vert(A) denote the set 
of vertices of A. For each vertex v G Vert(A), let 

/„ C {!,..., d} 

encode the set of facets that meet at u, so that 

i d Iv if and only if v E ai. 

Let ai^y, for i £ ly, be edge vectors emanating from v; concretely, assume that ai^y lies along the unique 
edge at v which is not contained in the facet cr^. (At the moment, the ai^y are only determined up to positive 
scalars.) In terms of the edge vectors, the tangent cone at a vertex v is 

Cy — {v + XjUj^y \ Xj > for all j}. 

The polar decomposition theorem relates the characteristic function of the polytope to the characteristic 
functions of convex polyhedral cones. As in the one dimensional case ()2.6f) . we cannot just consider the 
tangent cones, but we must make two modifications. First, we must "polarize" the tangent cones by flipping 
some of the edge vectors. Second, we must remove some facets. 

To carry this out, we choose a vector ^ gV* such that the pairings at^y) are all non-zero; we call it a 
"polarizing vector" and think of it as defining the "upward" direction in V. We "polarize" the edge vectors 
so that they all point "down" : for each vertex u of A and each edge vector a^.t, emanating from v, we define 
the corresponding polarized edge vector to be 

a.i^y if < 0, 



# 

i.v 



(3.1) a. 
Let 

ifv = [j e Iv I > 0} 

denote the set of "upward" edge vectors emanating from v, that is, those edge vectors that get flipped in 
the polarization process H3.1|l . The polarized tangent cone to A at w is obtained from the tangent cone Cy 
by flipping the jth edge and removing the jth facet for each j £ ipy-. 



(3.2) C*=lv+J2x,a*^ 

Recall that the characteristic function of a set A C M" is 

1 X £ A 



Xj > if j G Iy\ipy, and 

Xj > \i j £ Ifiy 



x(^ A. 



Theorem 3.2 (Lawrence, Varchenko). 

(3.3) 1a(x)=5](-1)'^"I1c*(x). 
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This decomposition was proved by Lawrence and Varchenko; see [O IGVI IV) . A version for non-simple 
polytopes appeared in |Haj . 

In preparation for our proof of Theorem 13 . 21 we introduce some notation. 

Let El,. . . ,En be all the different hyperspaces in V* that are perpendicular to edges of A under the 
pairing between V and V*. That is, 

(3.4) {E^\l<i<N} = {keraj- „ | v e Vert(A), j e /J. 

(For instance, if no two edges of A are parallel, then the number N of such hyperplanes is equal to the 
number of edges of A.) A vector ^ can be taken to be a "polarizing vector" if and only if it belongs to the 
complement 

(3.5) VI^V*\{EiU...UEn). 

The connected components of this complement are called chambers. The signs of the pairings ai^y) only 
depend on the chamber containing ^. 

Remark 3.3. For a Hamiltonian action of a torus T on a symplectic manifold M, one similarly obtains 
chambers in the Lie algebra t of T from the isotropy weights aj^p at the fixed points for the action. When 
Af is a toric variety corresponding to the polytope A, this gives the same notion of chambers as we have 
just described. For a flag manifold M = G/T, where G is a compact Lie group and T is a maximal torus, 
the isotropy weights aj^p are the roots of G, and the corresponding chambers are precisely the interiors of 
the Weyl chambers. 

Now suppose that ^ belongs to exactly one of the "walls" in (|3.4|) . Let e be an edge of A that is 
perpendicular to this wall. Let 

correspond to the facets whose intersection is e, so that 

i € le if and only if e C Ui . 

Let V be an endpoint of e. The edge vectors at v are aj^y, for j G /e, and an edge vector that lies along 
e, which we denote ae^y Note that 

{^,ae,v) = 0. 

Define 

(^e = {j e le I {taj,y) > 0} 

and, for each j G le, 

if {taj,y) < 
if {£„aj,y) > 0. 
Define the polarized tangent cone to A at the edge e to be 



a* = I 



Cf = { V + Xette + X^af, 



3 eh 



Xe e K, 

Xj > if i G /e\(pe, and 
Xj > if j £ ipe- 



Let v' be the other endpoint of e. We can normalize the edge vectors such that ae,t,' — —OLe,v and ctj.y—aj^yi G 
Mae,«. The set (p^ and the cone Cf' are independent of the choice of endpoint t; of e. 

Proof of Theorem \3./^\ Pick any polarizing vector ^ G V"^. Let v G Vert(A) be the vertex for which (^,t') is 
maximal. Then none of the a^^'s are flipped, and so Cj, = C.„. For any other vertex u G Vert(A), at least 
one of the OLj^yS is flipped, and so Cj^ n Cy — 0. So the polytope A is contained in the polarized tangent 
cone cf at v and is disjoint from the polarized tangent cone Cf for all other u G Vert(A), and the equation 
(|3.3|) . when evaluated at a; G A, reads 1 = 1. 

Suppose now that x ^ A. The set of vectors ^ which separate x from A, that is, such that 

(3.6) (C, x) > max ?/) , 

is open in V* . Choose a polarizing vector ^ G that satisfles H3.6() . Then x is not in the polarized tangent 
cone for any v G Vert(A). Equation (|3.3|l for the polarizing vector ^, when evaluated at x, reads = 0. 
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We finish by showing that, when the polarizing vector ^ crosses a single wall Ej in V*, the right hand 
side of (|3.3I) does not change. 

If Ej is not perpendicular to any of the edge vectors at the signs of (^,0;^,^) do not change, so the 
polarized tangent cone Cf does not change as ^ crosses the wall. The vertices whose contributions to the 
right hand side of change as ^ crosses Ej come in pairs, because each edge of A that is perpendicular 
to Ej has exactly two endpoints. 

For each such vertex v, denote by S^(x) and S^(x) its contributions to the right hand side of (|3.3|) before 
and after ^ crossed Ej. Let e be an edge perpendicular to Ej and v an endpoint of e. Let Sf{x) be the 
characteristic function of the polarized tangent cone Cf corresponding to the value of ^ as it crosses Ej . The 
difference Sv{x) — S^(a;) is plus/minus Sf{x). For the other endpoint v of e, the difference Sv{x) — S^(x) is 
minus/plus Sf{x). So the differences Sy{x) — S^(a;), for the two endpoints v of e, sum to zero. □ 

Remark 3.4. If we multiply both sides of H3.3|l by Lebesgue measure dx we obtain a formula for lAdx 
supported on A in terms of an alternating sum of the measures 1^^* {x)dx. This formula (which is a special 
case of "Filliman duality" |Fj) allows us to express the integral of any compactly supported continuous 
function / over the polytope in terms of its integrals over the cone . From the point of view of measure 
theory, the missing facets of Cf are irrelevant as they have measure zero; what is important is the change 
in the direction of some of the edges at a vertex and the sign (— l)!"^-! associated to the vertex. 

Remark 3.5. If the A is a polytope with non-singular fan, then (up to factors of 27r) the measure 1/^dx is 
the Duistermaat-Heckman measure of the associated toric manifold, and the vertices of A are the images of 
the fixed points of the torus acting on this manifold. In this case, the equality (|3.3|l multiplied by Lebesgue 
measure dx becomes a special case of the Guillemin-Lernian-Sternberg (G-L-S) formula. The G-L-S formula 
expresses the Duistermaat-Heckman measure for a Hamiltonian torus action in terms of an alternating sum 
of the Duistermaat-Heckman measure associated to the linearized action along the components of the fixed 
point set. See |(tLSI section 3.3]. This formula, in turn, can be deduced from the fact that a Hamiltonian 
torus action is cobordant (in an appropriate sense) to a union of its linearized actions along the components 
of its fixed point set. See |GGKI Chapter 4]. The role of the polarizing vector ^ G is played by a vector 
r] in the Lie algebra t such that the 77-component of the moment map $ is proper and bounded from 
below. 



4. Sums and integrals over simple polytopes 

Formulas for sums and integrals of exponential functions over simple polytopes appeared in jBVll IBV2| . 
In this section we deduce such formulas from the "polar decomposition" . 

We work in a vector space V with a lattice Vz- One may identify V with and Vz with Z", but we 
prefer to use notation that is independent of the choice of a basis. The dual lattice is 

(4.1) ^{ueV* \ {u, a) el for all a e Vj}. 

Remark 4.1. When our polytope is viewed as associated to a toric variety, the vector space V is the dual 
t* of the Lie algebra t of a torus T, the lattice Vz is the weight lattice in t*, and is the kernel of the 
exponential map t — > T. 

Let Cy be a simple convex polyhedral cone in V ^ that is, a set of the form 

n 

(4.2) C„ = {w + ^ XjUj \xj>Q for aU j} 

i=i 

where ai, . . . , a„ are a basis of V . Equivalently, we can write 

(4.3) Cy^{x\{u,,x)+\j>Q , i = l,...,n}, 

where ui, . . . , m„ are a basis of V* . We can pass from one description to another by setting ui, . . . , u„ to 
be the dual basis to ai, . . . , q;„ and vice versa. Geometrically, the vectors Ui G V* are the inward normal 
vectors to the facets of C^,, and they encode the slopes of the facets; the real numbers A; then determine the 
locations of the facets; the "edge vectors" ai generate the edges of C^,. 

A priori, the Mi's and the a^'s are only determined up to multiplication by positive scalars. To fix a 
particular normalization, we assume that the cone C^, is rational^ that is, that the a^'s can be chosen to 
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be elements of the lattice Vz, or, equivalently, the m^'s can be chosen to be elements of V^*. We choose the 
normal vectors Ui to be primitive lattice elements, which means that each Ui is in and is not equal to the 
product of an element of V£ by an integer greater than one. We choose ai, . . . , a„ to be the dual basis to 

Ui,. . . , Un- 

Although the u^'s are primitive lattice elements, they may generate a lattice that is coarser than V^*. The 
Q!i's then generate a lattice that is finer than Vz', in particular, the a^'s themselves might not be in Vz- See 
FigureEl The cone Ci, is non-singular if the UiS generate the lattice 1^, or, equivalently, the a^'s generate 
Vz. 



OL2 




Figure 5. 



Let dx denote Lebesgue measure on V , normalized so that the measure of V/Vz is one. Consider an 
exponential function f : V ^ given by 

(4.4) /(x) = e<«^^>, 

where ^ G is such that 

Re , ai) < for all i. 

Then the integral of / over the cone C^, and the sum of / over the lattice points in C^, both converge. 
If the cone Cy is non-singular and v is in the lattice Vz then the map 

{h, . . . ,t„) u + y^tj-gj 

sends the standard positive orthant 

M.l^{{ti,...,tn) >Ofor aUj} 

onto C„ and sends R" n Z" onto Cy CiVz- In particular, it takes the standard Lebesgue measure on R" to 
the measure dx on V. So 



00 /"OC 



f{x)dx = / • • ■ / e<«'''+S • • ■ dt„ 

c„ Jo Jo 

n poo 

7 = 1 "'O 



" 1 
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and 

oo oo 

c^nvk fei=o fe„=o 

n oo , 

(4.6) = e(«-) n E (^^^^"^0 ^' 

3 = 1 kj=0 
n 

A crucial ingredient in extending the Khovanskii-Pukhlikov formula to the case of simple polytopes is an 
extension of the formulas (|4.5() and (|4.6() to the case that the cone C„ is not non-singular. 
We associate to Cv the finite abehan group 

r = V^/span^{ui}. 

Note that the group F is trivial if and only if the cone Cy is non-singular. Also note that e^'^*^^'^^ is well 
defined whenever 7 G F and x G span^lai}. 

Remark. The toric variety associated to the cone is C"/F, where F acts on C" through the homomorphism 
F^ {S^y given by 7 (e^'^'^'^'^i^, . . . , e^'^^^'''""^). 

We have the following generalization of Formula (|4.6|l . Suppose that the vertex of satisfies v G 
span^jai}. Then 

n 

x^c^nVz ' ' 7er j=i 

Proof. The main step is to transform the left hand side of H4.7|l into a summation over elements of the finer 
lattice span2{Q;j}. For each x G span^jai}, 

(4.8) 7 ^ e^''''^^'^^ 

is a homomorphism from F to 5^, and it is trivial if and only if x ^ Vz,. Frobcnius's theorem asserts that, 
for a finite group F, the sum of the values of a non-trivial homomorphism F ^ S*^ is zero. It follows that, 
for X G span2{ai}, 

1 9^w^ ^\ I 1 if x G Vz; 



e 

F ^—^ otherwise. 



g27ri(7,x) 



(4.9) e^^'"^"''^ 
By (|4.9(l , the left hand side of (|4.7|l is equal to 

(4.10) E ^^'"^-mE 

2;GC„nspan2{Qj} 7er 

Writing 

X = V + kj ttj , 

this becomes 

n 00 
7Gr 3 = 1 k=0 

which is equal to the right hand side of (|4.7|l by the formula for the sum of a geometric series. □ 
We also have the following generalization of Formula (|4.5|l . 

(4.11) / e<«-)dx = e<«-).^n--l-. 
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Proof. We perform the change of variable x ^ v + X]j=i ^j'^j- Then x € Cy if and only if < = (ii, . . . , i„) 
belongs to the positive orthant M" . The inverse transformation is 

tj — {uj, X — v) ; 

its Jacobian is [V^ : span^Uj] — \T\. So 

Jc„ 1-1- 1 Js.1 

-1 ^ poo 

1 " 1 

□ 

To apply Formulas H4.7|l and (|4.11() to the "polar decomposition" of a polytope, we need to consider 
"polarized cones". Suppose that C ^ VJ? satisfies Re aj) ^ for all j. As in section |31 let 



for each j, let 

(4.12) 

and let 



fv = {j I Re{^,aj) > 0}; 



aj if j ^ ify 
~aj ifje^pv; 



Cl^\v + J2x 



> if j ^ ipv, and 

> if j e (^-u. 



Then the integral of / over Cj, and the sum of / over the lattice points in Cj, converge. We compute this 
integral and this sum: 



(4.13) 



/ e<«^->dx = e<«--> ~ ft V 



by 

by KT^ . 



Because v € span^jaj}, 



rrij G Z, 

Cj, n (span2;{aj}) = « + X! '^-J^i | ™j ^ if ^ V't', 

rrij > if j e (fly 



rrij G Z, 



TOj > for all j 



where 
(4.14) 

So 

(4.15) 
where 



c« n Vz - c" hift n 



c!,,shift = {vshift + X! I ^ for all j}. 
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We have 
(4.16) 

£cec{,nvz xGc" .hiftCys 



g{?, -"shift) 



1 1 

— V e2"<T'"=i""> TT = bv (im 



1 e~^'^*^'*''"j^e~^^^"j^ 



• 4^ y e'"^^'"^ TT , ■/ ^ ^ TT , , , by gH and Km 

76r i^¥'„ je<p^ 



(_l)lv.lg(«,f> . J_ ^ g2^»{7,f> JJ^ _ 



1 



(4.17) 



g27ri(7,Qj)g{5,Qj) 

■ ■ 7er 

1 

by applying the relation = — — to x ~ — 27ri (7. a,) — a,) for j £ ipy. 

1 — 1 — e ^ 

We can now reproduce Brion-Vergne's formulas for simple polytopes. Let A C be a simple polytope. 
Suppose that E satisfies Re {£,, aj v) ^ for all v £ Vert(A) and all j £ ly. With the notation of Section 

m 

■'^ •u6Vcrt(A) ■'^* 

i;GVert(A) ' ' jG/„ ' 

Similar formulas appeared in |BV1[ Proposition 3.10] and |BV2I p. 801, Theorem, part (ii)]. 
Also, 

xeAnVz t,eVort(A) aecfnVz 
a = y ^^^^^^ ■ 1^ y e'"^^'"^ TT . i \ r by Km 

(4.18) p 11 1 _ g27ri(7,Qj.„)g{4,aj,„) "« ' 

i)eVcrt(A) ' ' 7er„ je/„ 

i;GVert(A) ' ' jG/„ ' ^' ' 

where 

(4.19) Td.(5)= ^e^"<^--)TI i„g2..(7k.)e-^. ^"'^ ^ = {S,},ei.. 

7er„ je/„ 

Similar formulas appeared in |BV1I Proposition 3.9] and |BV21 p. 801, Theorem, part (iii)] 
Remark 4.2. We proved H4.17|l for ^ outside the real hyperplanes 

(4.20) Re(^,aj-^) , w £ Vert(A) , j £ ly 

in Vf^. However, the left hand side of H4.17|l is analytic for all ^ G V^, and the right hand side is analytic 
outside the complex hyperplanes 

i^'Otj.v) = 0. 

By analytic continuation H4.17|) continues to hold for all ^ outside these complex hyperplanes. Similarly, we 
proved (|4.18() for ^ outside the real hyperplanes ()4.20|l . but by analytic continuation is remains true for all 
^ outside the complex hyperplanes 

(4.21) (C, aj^y) = 2TTi {y, a,,y) , y£Vi. 

(Notice that these complex hyperplanes are contained in the real hyperplanes H4.20|l .'l We are grateful to A. 
Khovanskii for calling our attention to this approach. 
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Remark 4.3. Note that the function Td^(5') is analytic on the polydisk 

\ < bj , j e Iv 

in C^" , where 

bj = min { |27r (y, aj,^) | } . 

5. EULER MACLAURIN formulas FOR A SIMPLE POLYTOPE 

In this section we present Euler Maclaurin formulas for simple lattice polytopes in arbitrary dimensions. 
As in the previous section, we work with an n dimensional vector space V with a lattice Vz- Let A be a 
convex polytope in V with d facets, given by 

d 

(5.1) A = I {x,u,)+X,>0}. 

i=l 

The vectors Ui e V* are inward normal vectors to the facets, and they encode the slopes of the facets; the 
real numbers determine the location of the facets. As before, we assume that the slopes of the facets are 
rational, and wc choose the normal vectors Ui to be primitive elements of the dual lattice V^. We assume 
that the polytope A is simple, meaning that exactly n facets intersect at each vertex. We also assume that 
the Ai's are integers. 

Remark 5.1. If the vertices of A are lattice points, the A^'s are integers. If A has a non-singular fan and 
the Ai's are integers, then the vertices of A are lattice points. However, if A is simple but does not have a 
non-singular fan, its vertices might not be lattice points even if the A^ are all integers. 

For h near 0, the expanded polytope" 

d 

(5.2) A{h)^f]{xeV I {x,u,)+K + h,>0} 

is obtained from A by shifting the half-spaces defining its facets without changing their slopes. 

As before, we normalize Lebesgue measure on V so that a fundamental domain with respect to the lattice 
has measure one. The integral of a function / over the expanded polytope A(/i) is a function of hi, . . . , h^. 
For a polynomial P in d variables, the expression 

d _d_ 

dhi ' ' dhd, 

makes sense if the integral is a smooth function of h for h near 0. If P is a power series in d variables, the 
expression ()5.3|l makes sense if the resulting series converges. If P is an analytic function in d variables, we 
interpret the expression H5.3|) by expanding P into its Taylor series about the origin. 

In what follows, the function / can be taken to be the product of a polynomial function with an exponential 
function of the form e^^^^^ where ^ G is sufficiently small. 

Khovanskii-Pukhlikov's formula for a polytope with a non-singular fan. 

The formula of Khovanskii and Pukhlikov (see section 4 of |KP2p . translated to our notation, is the 
following formula. 



(5.3) P 



I 

h=0 JA{h) 



/i=A JA{h) 



where the polytope A is integral and has a non-singular fan. 

Finite groups associated to the faces of a simple rational polytope. 

The facets of the polytope A are 

CTi = {x e A I (wj,x)+Aj=0} , i^l,...,d. 

Because the polytope A is simple, each face _F of A can be uniquely described as an intersection of facets. 
We let Ip C {1, . . . ,d} denote the subset such that 
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The number of elements in Ip is equal to the codimension of F. The relation 

is order (inclusion) reversing. 

For each vertex v oi A, the vectors 

form a basis of V* . Let 

be the dual basis. The aj_„'s are edge vectors at v, that is, they point in the directions of the edges emanating 

from V. 

The vector space Vp normal to a face F is the quotient of V hy Tp = {r{x — y) \ x,y € F,r G M.}. Its 
dual is the subspace 

(5.5) Vp := spanjuj | j £ Ip} 

of V*. Let aj^p, j G Ip, be the basis of Vp that is dual to the basis uj, j € Ip, of Vp. 

To each face F of A we associate a finite abelian group Fi? in the following way. The lattice 

span^fuj I i G Ip} C VpCiV^ 

is a sublattice of n V^" of finite index. The finite abelian group associated to the face F is the quotient 

(5.6) F^ := (y; n V^) /span^lui | i G Ip}. 

For each 7 e Fjr and j e Ip, the pairing (7, aj^p) is well defined modulo 1, so 

g27ri(7,a;j,F> 

is well defined. 

Remark 5.2. The group Fir measures the singularity of the toric variety associated to A along the stratum 
corresponding to F. 

U F C E are faces of A, so that Ip ^ If, then {ui \ i € Jb} is a subset of {ui \ i G Ip}. Because these 
sets are bases of Vp and Vp , we have 

c y;, 

and 

Vp n span2{uj I i G Ip} = spanz{uj | i G Ip}. 
Hence, the natural map from F^ to F^ is one to one, and provides us with a natural inclusion map: 

a F CE then F^ C F^. 

We define 

(5.7) F^ = Ff\ U Tp. 

faces E such that EDF 

Finally, note that, for each face F and element 7 G Tp, the number e"^'^*^''''^^ is the same for all x G F; 
we denote this number 

g-27ri(7,F)_ 

Guillemin and Brion-Vergne formulas for a simple polytope. 

On any linear subspace AofV with rational slope we normalize Lcbesgue measure so that a fundamental 
domain with respect to the lattice AOVz has measure one. We shift this measure to any affine translate of 
A. Integration over each face F of A is defined with respect to these measures. 

For each face of A, let 

(5.8) F{h) = A{h) n f]{x\ {ui, x)+Xi + hi = 0} 

denote the corresponding face of the expanded polytope A{h). 

Guillemin gives an Eulcr Maclaurin formula for a polytope when the polytope is expressed as the set 
of solutions of the equation kiai + ... + kdCcd = M) ki,...,kd G M>o, for some fixed integral vectors 
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Qfrf, /i. (See Theorem 1.3 and formula (3.28) of [Gul.! When translated to our setup, his formula 
becomes the following formula. 



f5 9) V f = V — V e^"^'^'-^^ TT ^ TT - 

AnVz F ' ^' ^gr^, j^iF jeiF 

where the polytope A is simple and is given by (|5.1I) where all the A^'s are integers. 
Finally, the formula of Brion and Vergne in our notation is 



/i=0 



F{h} 



a 

^27ri(7,F> 77 Y\ ^ 

unvz F ^grj, j^ip " ~ jell 



(5-10) E / = E E ^'"'^^'^ n n r 



A(/0 



h=0 



See jBVll Theorem 2.15] (where A is simple and integral), and, more generally, |BV2| . 

Remark 5.3. If the polytope A is integral (which is a stronger requirement than the assumption that the 
Ai's be integers) then e^'^*^''"'^^ = 1 for each face F and each 7 G Tp. 

We now give a self-contained statement and an elementary proof of the Guillemin-Brion- Vergne formulas. 

Theorem 1. Let V be a vector space with a lattice Vz- Let C V* be the dual lattice. Let A <Z V be a 
simple rational polytope with d facets. Let ui, . . . ,Ud G be the primitive inward normals to the facets of 
A. Let Xi, . . . , Xd be the real numbers so that 

d 

A= f]{x I {u„x) + X, >0}. 
1=1 

Suppose that Xi, . . . , Xd are integers. For h — (hi, . . . , hd), let 

d 

Mh) = I +X, + h,> 0}. 

1=1 

For each face F of A, let Ip C {1, . . . , d} be the subset such that F consists of those x G A for which 
{ui,x) + Ai = for all i G Lp. Let ai^p, for i G Ip, be the basis of V/Tp = y/M(_F' — F) that is dual to 
the basis Ui, i G Ip, of Vp = Tp. (In particular, if v E Vcrt(A) and i E ly then Ui^y are the edge vectors 
emanating from v.) Let 

Tp = (y; n yz*)/spanz{ui | i G /f} 
be the finite group associated to the face F, and let T"p — T p'\\_^T p, where the union is over all faces E 
such that E D F. Let 

Td.(5„...,^.)=E E ^'-^-'^'^ n n i_,2..(l.),-s. - 

F 7er5, ]<^If jeip 

Let f : V C be a quasi-polynomial function, that is, a linear combination of functions of the form 

fix) =p(x)e<«'-> 

where p: V C are polynomial functions and where the exponents G satisfy 

(5.11) I (e,aj-„) I < 27r| (y,aj-„) I 

for each vertex v, each edge vector aj^y, j G /„, and each y G such that {y,aj y) ^ 0. (Each of the sets 
{{y,aj y) I y G } discrete, so the set of^'s that satisfy (b) is a neighborhood of the origin in V^.) Then 

(5.12) ^ /(.,.Td.(^....,|^: 



f{x)dx. 

h=a JA{h) 



Remark 5.4. The right hand side of H5.12(l is an infinite sequence. The theorem asserts that this sequence 
converges to the left hand side. In Appendix^we show that this convergence is uniform on compact subsets 
of fTTTIl . 

The proof of Theorem uses the following characterization of F^. 
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Lemma 5.5. Let F be a face of A. 

(1) If ^ £ Tp and v is a vertex of A such that v £ F, then 

g2.z(^.o,.„) ^ fe'"*<^'"--> foralljelp 

1 1 for all j G ly-^Ip. 

(2) Ifj e If and 7 G T^, then eS'^^^T."^.^-) ^ 1. 

(3) For 7 G r^., 

7 G if and only if e^'''^^'"^-^^ 7^ 1 /or all j £ Ip- 

Proof Let G n V£ be a representative of 7 (see (ESJ). Then, by definition, e^'^^^'^'^^vp) = e^'^*^!''"^.^-) . 
Because y G V^, and by (|5.5|l . there exist real numbers Oj, for j G /f, such that y — J2j£ip Then 

(5.13) {y, cxj^p) = o,j for aU j G Ip. 
Defining aj = for j G Iy\Ip, we also have y = Sje/ '^i'^j' ^'^'^ 

(5.14) {y,aj^v) = aj for all G Iv 
In particular, 

(5.15) {y, aj.v) — for all j G Iv\If- 

Part (1) follows from fO^ . (^33, and 

Fix j G If. Suppose e^^*^''''"^-^^ = 1. Then we can choose a representative y G fl of 7 such that 
{y,aj,p) — 0. Writing y — "^i^j^ aiui, we have aj = {y,aj,p) — 0. Let E be the face of A such that 
Ie = Ip^ij}- Then y = X^ie/^ azit;, so, by H5.5(l . y is in FJ, and so 7 G Fs- In particular, by H5.7|l . 7 ^ F^. 
This proves Part (2). 

Let 7 G Ff. By Part (2), if 7 G F"^ then e^'^'^T^"^.^) ^ 1 for all j G /f- Conversely, suppose that F ^ F^. 
Then, by (|5.7|l . there exists a face £^ such that 7 G Ff and E D F. Let j G Ip-^Ip. Let w be any vertex of 
F (and hence of E). Then e^'^*^''''"^'*'^ = e^''*^''''"^'''^ = 1, where the first equality follows from Part (1) for 
the face F, and where the second equality follows from Part (1) for the face E. This proves Part (3). □ 

Claim. For each v G Vert (A), 

(5.16) TdA(S'i, . . . , Sd) = Td„({5j}je/J + muhiples of for j ^ /„. 
Proof. Recall that 

5, TT ^ 



Td.(5„...,5.)=E E e"^^^'"^ n n T 



_ g2iri{7,aj,F)g-Si ' 



By Part (2) of Lemma [5. 51 for each 7 G T^p and j Ip, 



S 

— a multiple of Sj 



Td.(5„...,5.)=^ E, E^""^'"\ n, T^ni 

By 



j ^ /i, , and becai 

5*^ -p-|- Sj multiples of Sj 



Because v ^ F implies that there exists j G Ip such that j ^ Iv, and because j — rr^- = 1 + a multiple of Sj , 



F such that ^^F^ j^^v^Ir j^Ip 

«GF ^ 



F„ = □ F^. 



F such that 
■U6F 



Also, e2^*<'>''-^> = e^'^'^'^^''^ whenever v G F. By this and Part (1) of Lemma ICT 

TdA(5i, ...,Sa)^Y. ^"^^^""^ n i„g2..(^,„)g-5, + ^^iltiple^ °f fo'- ^ ^ 
By the definition (|4.19|) of Td^,, this exatly shows (|5.16|) . □ 
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Proof of Theorem n Khovanskii-Pukhlikov approach: 

Let C M'^ be the set of all h € M.'^ that are sufficiently small so that the polytope A{h) has the same 
combinatorics as the polytope A (i.e., the same subsets Ip C {!,... ,d} correspond to faces). The vertices 
of the expanded polytope A(/i) are then 

(5.17) v{h) — V — hjUj^y, 
and we have 

(5.18) ^e<«-W)=-(e,a,,.)e<«-W> 
for any j iz Iv Let 

I^{h,v)^ f e<«'^>da; and 5a (C) - e<«'^>. 

By BTTIl . 

(5-19) ^A(^e)= E ^''■''-'■^^W-jrhi 

^6Vcrt(A) I "Ije/^ 

for all e and ^ £ that lies outside the complex hyperplanes 

(5.20) (^,aj-,}=0 , uGVcrt(A) , j G 
By (jTTH|l . 

(5.21) 5a (0 = E e<«->Td„({- (C, a,v«)},g/J ' pf^ I] - /F^' 

^eVert(A) ' "'iG/^ 

for all f G that he outside the hyperplanes (|4.21() . See Remark 
By (|^T7|l . and (EHH), 



(5 



22) TdA(^,...,^)lA(/.,C)= E Td„((Aj )e<«-(")> . n -TF^ 



So the partial sums of the series H5.22|) are 

(5.23) ^p„^^({_(^,«^.„)}^.^,Je(«.''W). 1 [] -TF^ 

where Pm,v are the Taylor polynomials of Tdt,. By Remark lOl Pm,v{{Sj}j,^i^) converges to Td„({5j}je/^) 
uniformly on compact subsets of the polydisk 

(5.24) {S G C^^ I \Sj \ < bj,y for aU j G Iv}, 
where 

(5.25) bj,v — mill 27r (y, aj.u) . 

(y,a3,,)7^0 

Because the functions 

are (continuous, hence) bounded on compact subsets of the set of ^'s that lie outside the complex hyperplanes 
given by H5.20|l . the partial sums 15.23(1 of the scries 1(5.22(1 converge to 

E Td.({- {^.a,.)}^,Je(^-(^^^ ■ -L H 

uniformly on compact subsets of the set of (ft,, ^) G f2 x such that is outside the hyperplanes ((5.2U() and 
satisfies 

(5.26) I aj^y) I < bj^y for aU v G Vert(A) and j G /„. 
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Setting = 0, by H5.21|l . we get 



h=0 



and that the left hand side of (|5.27|) converges to the right hand side of (|5.27|) uniformly in ^ on compact 
subsets of the set of ^ G T^; that lie outside the hyperplanes H5.20|l and in the set (|5.26() . 

However, the right hand side of H5.27|l and the partial sums of the left hand side of H5.27|l are analytic in 
^ for all ^. 

Recall that, as a consequence of Cauchy's integral formula, if g^iO is a sequence of complex analytic 
functions on an open subset U of C", g{£,) is a complex analytic function on gv(£,) converges to g{S,) in 
U\E where _E is a locally finite union of complex hyperplanes, and this convergence is uniform on compact 
subsets of U\E, then g,y{^) converges to g{^) for all £, E U, uniformly on compact subsets of U. 

It follows that (|5.27|) holds for all ^ G V,^ that satisfy (|5.2t)|) . and, moreover, the left hand side of (|5.27f) 
converges to the right hand side uniformly in ^ on compact subsets of (|5.26f) . This gives the Euler-Maclaurin 
formula for exponential functions e^^'^^ for all S, in the set (|5.26|l . It also shows that the limit on the left 
hand side of 15.27|) commutes with differentiations with respect to ^. Applying such differentiations to the 
left and right sides of (|5.27|) . we get 



whenever P{x) is a polynomial and ^ is in the set (|5.26|l . 

In particular, for ^ = 0, we get the Euler-Maclaurin formula for polynomials. 

Proof of Theorem ^ for polynomial functions Brion- Vergne approach: 

The terms in H5.19|l and (|5.21|l arc functions of ^ whose products with Yij v extend to analytic 

functions near ^ = 0. Comparing the Taylor expansions in ^ of the left and right hand terms of these 
products, we get 

k+n 



(5 28) f^^dx- V iii^.^n 

^ ' Ja k\ ^ (fc + n)! |r„| -11 

•'^ ueVcrt(A) ^ ' I "I jeh, 



and 



("8) E E (=«->™.(i-«.»„)wJ)"""^p^n 

where the superscript {k + n) denotes the homogeneous term of degree fc + n in ^. 
Recall that the vertices of A(/i) are 

(5, v{h))^ is a polynomial of degree k in the /ij's that only depends on hj for j G I^. For all j e /„, 



dhj kl (fc-1)! 
So for any homogeneous polynomial T(-) of degree £ in the variables Sj, j E 1^, 



(5.30) T ( {^},,,. 



h=0 



= ^ ) (fc^£)! ■ 
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We have 



d d \ f {^,x)' 




«eVert(A)0<£<fe+n / ^K + n «j. I^^l^-g^^ ^i;,"^,^ 



fc! ■ 



When h — 0,hy H5.29|l . this is equal to 

E 

6. The Stokes formula for polytopes and the Cappell-Shaneson algebra 

Khovanskii and Pukhhkov work with derivatives associated to the "expansion" A{h) C of the 
polytope. In this section we give two results that relate such derivatives to differential operators on V. The 
first result, Proposition 16.11 is the Stokes formula for polynomials. The second result, Proposition 16.21 is 
that integration over faces can be replaced by differentiations with respect to corresponding hj^s. A similar 
argument appears in |BV2[ Section 3.6]. We use these results to define the "Cappell-Shaneson algebra", a 
formalism used by Cappell and Shaneson to express their formulas. These results play a key role in relating 
the Cappell-Shaneson formula to the original Khovanskii-Pukhlikov formula; we do this in Section [7| The 
two results can also be used to derive the Euler Maclaurin formula for polynomials from the formula for 
exponentials, as we do in Section IbI 

The Stokes formula for polytopes. 

Proposition 6.1. Let V be a vector space with a lattice Vz- Normalize Lebesgue measure on V so that the 
measure of a fundamental domain for the lattice Vj, is one. Let A be a rational polytope in V . Let ui, . . . ,Ud 
denote the inward normals to its facets, normalized so that they are primitive elements of the dual lattice V£ . 
For any v gV, let Dy denote the directional derivative in the direction of v. Then, for any f G C°°{V), 



(6.1) 



Proof. The formula is an immediate consequence of the general Stokes formula. 
Alternatively, it follows directly from 



A{hi + {ui,v} ,...,hn + (u„, v}) = A{h) - v. 

□ 



Integration over faces 



Proposition 6.2. Let A be a simple polytope and let F be a face of A. Let A{h) be the expanded polytope 
and F{h) the corresponding face of A{h). (See (|5.1|) . (|5.2|) . and (|5.8|l . ) Then, for any smooth function 
f G C°°{V), the integral of f on A{h) is a smooth function of h for h near 0, and 

(6.2) / f=\^F\\[^! f 

where Tp is the finite abelian group associated to the face F. (See (|5.6() . ) 
In particular, 

(6.3) J^...^ f f^o ifa.,n...na,, =0, 
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and 

w [ ^= I ^- 

Ori-i JA{h) Jai{h) 

Proof. Choose a polarizing vector ^ G such that if u £ Vert(A), v ^ F, and x E F, then x) > v). 
(For instance, we may take ^' such that the restriction of the hnear functional (^', •) to A attains its maximum 
along the face F, and take ^ to be a perturbation of which is in V^.) Then the edge vectors that are 
based at a vertex of F but are not contained in F (that is, aj^y for v G F and j £ Ip) are not flipped in the 
polarization process H3.1|l . 

Let Pp denote the affine plane generated by the face F. After possibly multiplying / by a cut-off function 
which is equal to one near F, we may assume that {$,,x) > {^,v) for every vertex v which is not in F and 
every x G Pp D supp(/), where supp(/) is the support of /. 

Then for every vertex v which is not contained in F we have 



f = o. 

Similarly, 

(6.5) / / = if V ^ F and h is sufficiently small, 

iPf.(h)nc*(/i) ■ 

where Pp(h) is the affine plane generated by F{h) and where Cf{h) are the cones that occur in the polar 
decomposition 

(6-6) l^(,)^^(-l)l'^"ll^*(,). 

V 

From (|6.5|1 and (|6.6|1 we get 

JF{h) JPrih) JPF{h) JPF^nCtih) 

It remains to show that 

(6.7) / / = ir^^i U^r [ f 

^p^,(^)nc*(/i) ^^j^ on-t Jc*(h) 

for each v £ Vert(A) such that v E F. Assume without loss of generality that hj = for all j G Iy\Ip. 
Let 

Tp = {r{x - y) \ x,y E F} 
denote the tangent space to the face F. Consider the affine change of variable map 

ip: Tp xR^'' ^ V 

given by 

(p{y,t) =v + y+Y^ tjaj.y. 

Let 

F*^ J2 R+«t and R[-{h)= Hi-h,,^). 

The map if sends F* X M^^ onto Cfih) and sends F* x {{~hj)jeiA onto F*{h). Lebes gue measure m 
F^ C Tp is normalized so that the measure of a fundamental chamber for the lattice Tp n Vz is one. Clearly, 

/ / My, ~h)) dy^Yl^ f f (^(y, t)) dydt. 

To conclude 1)6.71) it remains to show that |dct dip\ ~ jpj^- 
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The map dip sends the subspace TF x {0} C TF x M^^ to the subspace TF C V and respects the 
lattices in these subspaces. So its determinant is equal to that of the induced map on quotients. Recall that 
V/TF — Vp- The induced map on quotients is the map 

Tp: M.^" Vf 
given by ^((tj)je/F) = Eje/j. ^j^h^- inverse, 

V': Vf 

is 

V;(/?) = ((u„/3)),e/.. 

The dual ip* : MJ'' — > Vp sends the standard basis element Cj to Uj for each j E Ip- Finally, 

detd.^ = det^= (det-0)"^ = (det-^*)"^ = [spa,n.^{uj} : V^nF^*]"^ = iTi^P^ 
as desired. □ 

By H6.2|l . the formulas of Guillemin 1)5. 9|l and of Brion-Vergne (|5.1()|l are equivalent. 

The Cappell-Shaneson algebra 

Let F be a vector space with a lattice Vz and A C V a simple lattice polytope. Let V denote the ring 
of infinite order constant coefficient differential operators on V. Consider the algebra I?[[[cri], . . . , [ad]]] of 
power series in the formal variables [tXi], corresponding to the facets, and with coefficients in V. A general 
element of this algebra can be written as 

d 

a i—1 

where Pa E T> for each a: {1, . . . ,d} Z>o. Each such element A defines a linear functional / A which 
associates to each polynomial f on V the number 



(6.8) 



pM)- 



J A{f) is also defined for every smooth function / for which the right hand side of IjG.Sfl is absolutely 
convergent. By H6.3|l . 



if (Ti^ n ■ ■ ■ n (Ti^ 



(6.9) JW^^]■■■['J^.] 

By (EH) and 

/ + E [^'] j = ^ V eV. 



(6.10) 

Following Cappell and Shaneson, we consider the quotient Q{A) of 2?[[[(7i], . . . , [(T^]]] by the equivalence 
relations 

(6.11) K]---[a,J=0 

for each set (7^^ , . . . , at^ of facets where fl . . . fl Cij. = 0, and, for each v E V, 

d 

(6.12) D,+J2{v,u,)[<j,]=0. 

1=1 

We call (5(A) the Cappell-Shaneson algebra. (Cappell and Shaneson denote this algebra Q(S), where S is 
the corresponding fan.) 

Remark 6.3. The quotient of the polynomial algebra 2?[[cri], . . . , [ad]] by the relation (|6.11|) is the Stanley- 
Reisner ring (the face ring) of A with coefficients in V. 
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By and the definition of / A{f) descends to tfie Cappell-Sfianeson algebra (5(A): for T £ Q{A) 

that is represented by A S I?[[[(Ti], . . . , [crd]]], we can define 



(6.13) J T{f) = J A{f). 



Remark 6.4. Cappell and Slianeson do not consider variations A{h) of the polytope A. They consider the 
free D-module P{A) with basis elements [F] corresponding to the faces F of A and the D-module map 
p: P(A) Q(A) defined by 

pm)^\rF\l[[a.]. 

For T = p{n), with = 6 P{A), Equation implies 

(6.14) J^(f) = T.J/Pf- 

Cappell and Shaneson define J T{f) by H6.14() . Comparing with 16.1 51 we see that this is well defined. 
Define the summation functional by 

sif) = E /■ 

AnVz 

Using the Cappell- Shaneson algebra, the Khovanski-Pukhlikov formula (|5.4|l for an integral polytope with 
non-singular fan reads 



S 



d r 



i=l 



and the Guillemin-Brion-Vergne formula for an integral simple polytope reads 



o2iri{7,Q!j_F>g-[CTj] ' 

1. The Cappell-Shaneson formula 

The main differences between the formulas of Cappell and Shaneson and those of Khovanskii-Pukhlikov, 
Guillemin, and Brion-Vergne, are these: 

• The latter authors work with expansions of the polytope. Cappell and Shaneson work with derivatives 
of the function on V . 

• Cappell and Shaneson express their formula in terms of what we call the "Cappell-Shaneson algebra" . 

• Cappell and Shaneson derive their formula for the sum X^AnVz ^ f^'om formulas for the weighted sum 

/ -J \ codimi^ 
AnVz _F ^ ^ rcl-int(F)nyz 

Cappell and Shaneson's exact formulas for simple lattice polytopes, when applied to polytopes with non- 
singular fans, become the following formulas: 

Theorem 2. Let A fee an integral lattice polytope in a vector space V with a lattice Vz- Let f he a poylnomial 
function on V . Then 



(7.1) E / / ^(/) 



AnVz 



(7.2) Yl nf), 

interior(A)nVz 

and 

(7.3) E / - J E /= l\^T + f){fl 



2 •' / 2 

AnVz SAnVz 
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where 

r = V TT M TT [^']/^ 

4^ il 2 ii tanh([ai]/2) 

anrf 

T = 'S^C-n™'^™-'^ TT 1^ TT I'^']/^ 

Z^^ ^) ii 2 J-i tanh([(7i]/2) ■ 

The relation of these formula to the Khovanskii-Pukhlikov formula goes through similar formulas that 
apply to a polytope with some facets removed. We first need a corresponding polar decomposition. 

As before, let A be a polytope with d facets. For a subset L C {1, . . . , d}, denote by the set obtained 
by removing from A the facets (Tj, i £ L. In particular, for L = 0, A^ = A, and for L = {!,..., d}, 
A^ = interior(A). 

Recall that A = Hi D . . . (1 where each Hj is a half-spaces whose boundary is the afBne span of the 
facet aj. Let 

I interior (f/j) ifj^L. 



Then we have 



3 

Fix a polarizing vector ^ for A. Recall that this determines a subset of for each vertex v. Consider 
the cones 

where 

„#,!, ^ / Hi if 3 G Iv\V>v 
^> \{Hj^y ifje^.. 

We have the following polar decomposition for A^: 
Proposition 7.1 (Polar decomposition with some facets removed). 
(7.4) l^.(x)=^(-l)l'^''ll^*,.(x). 

V 

Proof. We shift the bounding hyperplanes of Hj inward or outward according to whether j € L or j ^ L. 
That is, we shift by an h which belongs to the set 

Orth^ := {{hi, ...,hd)\hj<0 for j e L and hj > for j ^ L.} 

We have the pointwise limits 

1^l{x) = lini lA(ft) (.'£■) 

heOrth^ 

and 

1^^.l{x)= lim lc#(ft)(a;) 

/leOrth^ 

for all X. The proposition follows immediately from these limits and from the polar decomposition theorem 
for A{h). □ 

We have the following variant of the Khovanskii-Pukhlikov formula for a polytope with some facets 
removed: 

Proposition 7.2. Let V be a vector space with a lattice Vz- Let A be a lattice polytope in V with facets 
CTi, . . . , (Td and f a polynomial function on V. Let L c {1, . . . ,d} be any subset. Suppose that A is a polytope 
with a non-singular fan. Then 



E/^/nS^;nT^(/). 



EXACT EULER MACLAURIN FORMULAS FOR SIMPLE LATTICE POLYTOPES 



25 



Proof. The proof follows exactly the same lines as the proof of the Khovanskii Pukhlikov formula, (|5.4|l . 
usmg the polar decomposition (|7.4|) for A-^. We leave the details to the reader. □ 

We derive the following formula for the weighted sum: 

Proposition 7.3. Suppose that A is an integral polytope with a non-singular fan and f is a polynomial. 
Then 

/n J?^</) - 

Proof. Consider the left hand side of equation H7.5|l : 

(7.6) ffjJ^iI^m. 

^ > i litanhK-/2]^^^ 

Since 



tanh(i:)/2) ' ' 'e^f^-e-D/^ 2' ' 1 - g-^ ' 

()7.6|l is equal to 

Expanding, this becomes 

2) 51 / n j !_^e-k,] n I Jg^-k.] 

' ^ LC{l,...,d}"' jeL j^L 

which, by Proposition l7.2l is equal to 

For each face F of A, the relative interior of F is contained in A^ if and only if /j? n L = 0. The number of 
subsets L which satisfy this condition is 2'^"'=°'*™^. Therefore, H7.7|l is equal to 

/ -| \ codimi^ 
F ^ interior(F)nVz 

which is the right hand side of (|7.5|) . □ 
Lemma tl.4\ relates weighted sums to non- weighted sums. For a face F of the simple polytope A, let 

E'/ 

AnVz 

denote the weighted sum with respect to the affine span of F, that is, 

/ \ dim _F— dim i? 

E'/=eQ) e /• 

AnVz E ^ rel-int(_E)nyz 

E<ZF ^ ' 

Here, the faces of F are exactly the faces i? of A that are contained in F, and the exponent dimF — dimi? 
is the codimension of E in the affine span of F . 

Lemma 7.4. 

codimi^ 



and 

codimi^ 



Ez-EG E'/, 

/ ^ \ codimi^ 

E /-eR) E'/ 



X 2, 

interior(A)nVz -F ^ -FnVz 
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Proof. 



\ codimi^ / 2^ \ codimi^ / \ dim i^— dim i? 



±7 T.'f-T.^ 



V 2 / ^ ^ \ 2 1 ^ V 2 / ^ 

F ^ / FnVz F ^ ^ E ^ rel-int(B)nyz 

_ECF 



/ -| \ codimF / 

(7.8) ^(±l)codimF ^ ^ 



F 
,FDB 



rol-int(F)nyz 



Because A is simple, 



1 = 2™^™^ and (-i; 

F F 
FDE FDE 



codimF' 



E CA 

1 E = A. 



Substituting this in (|7.8|l gives the lemma. □ 
We are now ready to derive the Cappell-Shancson formula. 

For each face F of A, we have z G /i? if and only if ai D F. For i ^ Ip, the intersection ai H F is either 
empty or is equal to a face of A which is a facet of F. We denote 



Since 



and by (j6.9|l . we get 



Sf = {« I n F is a facet of F}. 



tanh([a-i]/2' 



= 1 + a multiple of [ai], 



(7.9) ^=En¥n 



F .,1. 2 tanh([a.]/2) 



and 
(7.10) 



^^^(_,)codimF-QM -Q 



2 J- J- tanh([(Tj/2) 



as elements of the Cappell-Shaneson algebra Q{A). By (|7.9|l and H6.2|l . followed by Proposition 17. 31 applied 
to the face F, and further followed by Lemma 17.41 we get 



r / 1 \ codimF p r 1 /O / 1 \ codimF 

Similarly, from HT.lOfl we get 

F ^ ^ "'-^ jGSf U ij/ ; p \ / pf^y^ intGrior(A)nV2; 

This proves H7.1|l and H7.2(l . The equality (|7.3|l clearly follows from these. 

Remark 7.5. We expect that a similar argument will show that the Cappell-Shaneson formula for simple 
polytopes is equivalent to the Guillcmin-Brion-Vergne formula. 

Remark 7.6. Formulas with more general weightings have been developed in |AW| . 



odimF /. r 1 /r» / -i \ codimF 
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Appendix A. Relation to remainder formulas 



Another exact formula for polynomial functions on simple polytopes appeared in our recent paper |KS W2| . 
There we proved an Euler Maclaurin formula with remainder for simple polytopes and gave estimates on 
the remainder. From this we deduced an exact formula for polynomials directly, without passing through 
formulas for exponential functions. Let us describe our exact formula from |KSW2| in our current notation. 

Let A be a complex root of unity, say 

A^ = L 

Define a sequence of functions Qm,\{x) on E recursively, as follows. For m — 1, set 

A 



Qi,x{x) = YTT^X! ^"■^["^"+l)(^)• 



Given the function Qm-i,\{x), define the function Qm.\{x) by the conditions 

dx 

Consider the polynomial 



.\{x) = Qra-i,x{x) and 



M*'-^(S') 



1 A 

t: + 



\{x)dx = 0. 



2 1 - A, 

Let V he a, vector space with a lattice Vz- Let 

A = {x \ (ui, x) + Hi>0 , i = l,...,d} 

be a simple lattice polytope in V, where ui, . . . ,Ud £ V* are the normals to the facets of A, normalized so 
that they are primitive elements of the lattice V£. Let 

A(/i) — {x I {ui, x) + + hi > , i = 1, . . . , d}. 

For a face F of A, an element 7 of Fi^, and an index I < j < d, let 



A, 



1 



and consider the differential operators 



j e If 



""dhi 



Let 



AnVz F rcl-int(_F)nVz 

summing over the faces of A. Then for any polynomial function / on A, for sufficiently large k, 



(A.l) 



k 

■'r,F 



f 



A{h) 



h=0 



For comparison, the Euler Maclaurin formula (|5.12|l for simple lattice polytopes can be written as 

E/-EET../ / 

AnVz F T,er5, ^^''^ 

with 



h=0 



-1,F 



_d_ 

dhi 



and 



S 



1-Ae-S' 



A similar argument (see below) gives 
(A.2) 



E7 = EE 

AnVz F ^grj. 



A(/i) 



h=0 
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ith 

d 



and 



As observed by Michele Vergne |Ve2| , the equivalence of Formulas (|A.1(I and l)A.2p is seen from 
Lemma A.l. M'^''*'(S') is the kth Taylor polynomial ofL'^{S). 

We complete this section by giving the proofs of (|A.2p and of Lemma FA. II 

Proof of (|A.2|I . We have the following analogue of 1)4.7(1 : 

xec^.nVz I I 7er j=i 



Indeed, applying 



1 1 + Ae ^ 1 ^ _o ,2-2'? 
- ^ Ae + A^e 



2 1 - Ae-^^ 2 

to A = e27ri{7,ti) g-s _ g(C:Qj> g^jjfj rearranging the terms, the right hand side of (|A.3|I is equal to 

(fci,...,fe„)ez|o 76r 

which, by 1)4. 9|l . is equal to the right hand side of 1A.3|I . 

From this we get the following analogue of 14.16|l : 
(A.4) 

1 "ill 27ri(7,Q".) ii-a") 

1 'sr^ 9^,7^ „^ T-r 1 1 + e^'^^^T'^^^e^^'"^^ -rT 1 1 + e-2"(7,",)e-(?'".) 



|r| .J- J- 2 1 - e2"(7,a,>e{?,a,> 11 2 1 - e-2"(7.",)e-{?,".> ^ 



II 'I _ ^ ^ - ^ - 

7er jev 

1 V- 9.,/. „\ TT 1 1 + e2^*<T,a.)e(«.".) 



g27ri(7,Qj)g(^,Qj) 

■ ■ 7er 

1 + 1 + e"^ 

by applying the relation — — to x = — 27ri (7, a,) — a,) for j d tp^. 

1 — 1 — e ^ 

Let IX (a^) denote the weighted characteristic function, given by IX (2;) = (^1^'^°'^""^ jf 2; ijgs in the relative 
interior of a face -F, and IX (a^) = if a; ^ A. Define 1^(2;) in a similar manner whenever C is a convex 
polyhedral cone. We have the following analogue of 1)3.3(1 : 

(A.5) lX(a;)=^(-l)l^"ll-,(:.). 

V 

This can be proved directly (see |KSW2[ section 3]), or it can deduced from ((7.4(1 using the formulas 
^a{x) = ^ 1a4x) and 1^, (a;) = -^l^j,!- (x). 

LC{l,...,<i} 
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From this wc get the following analogue of H4.18|l : 
x&AnVz ■ueVort(A) xe:C*nVz, 



(A.6) 



g27ri(7,aj.„)g({,aj,„> 



by (ESI 
by ^EM 



uGVcrt(A) 

where 
(A.7) 

Note that Lt,(S') is analytic on the polydisk {\Sj \ < bj, j G /„} that is described in Remark l|4.3|l . 
The operator that appears in ljA.2p can be written as LA(gf^, • ■ • , gf^) where 

d 

We have the following analogue of H5.16|l : for each v G Vcrt(A), 

(A.8) La (5i , . . . , Sd) = L„ {{S, ) + multiples of S, for j ^ I,. 

This is shown exactly like H5.16|l . using the facts that if A 7^ 1 then L^(S') is a multiple of 5 and if A = 1 
then L-^(5) = 1+ a multiple of S. 
By (|5T7jl and (|5l^ . 



(A.9) 

u6Vert(A) ' ^ ' 

"2)1 follows from (|A.9p and (|A.8|) by the same arguments as in the proof of Theorem ^ in section [S] □ 
Proof of Lemma \A.l[ Suppose that A 7^ 1 and A^ = 1. An argument similar to those in section [21 gives 



d 



dhi 



d 



hi=h2=0 J ~hi 



N + h2 



f{x)dx 



dh2, ^ 

V(0) + A/(l) + A2/(2) + . . . + A^-V(A^ - 2) + X^-'fiN - 1) + l-f{N) 



for all polynomial functions f{x). Indeed, direct computation of ^ + Ae^ + A^e^^ + . . . + A^^^e'^^^^^ + i, 
followed by multiplication by ^ and taking the A'^-th degree term in the Taylor expansion, gives the following 
analogue of H2.12|l : 

(A.ll) e Y '(e)-(L^"(0-e«^-L\-0-l)^''^'^ 

xe[o,N]n7, 

whenever Ae^ ^ 1, where the superscript {N + 1) denotes the N + 1th term in the Taylor expansion. On the 
other hand, 



/i=0 ''-h 



N+hi (c^\N 



m 



-dx 



d 



hx=Q 



(iV + l)! ^dh2 



N+l 



H2=o (N + iy. 



by (ITTTTl 

{N + l) 



by (E3 



By direct computation, the first Taylor coefficient of L'*'(S') is | + and that of L'^ (5) is 5 + 



i — yta- Let am denote the mth Taylor coefficient of L'*'(5'). Since L'*'(— S") = L"** (5), the mth Taylor 
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coefRcient of ^ (5) is (— l)'"am- Taking F{x) to be a polynomial of degree < fc + 1 and f{x) — F'{x), the 
left hand side of (|A.10p becomes 

On the other hand, the right hand side of (|A.10p is equal to 

1 A 11 A 

(^•13) (9 + T-A + ( 9 - T-A + T-A E + 1) - /(")) 



n=0 



1 A 11 , ^ Af /"^ 

2 i-A 2 i-A -'0 

for any fc > 2, by repeated integration by parts as in the proof of Proposition 27 of | KSW2 |. Recalling that 
Qm.xiO) — Qm,\{N) and that Qk ^\{x)dx = 0, taking f{x) — F'{x) where F is polynomial of degree 
< fc + 1, the right hand side of (|A.13p becomes 

+ 1^)^^°^ + T^)^(^^ + ^ (-i)"Q™.a(o) (f(")(7V) - ^^(™)(o)) . 

T?X = 2 

Comparing this with ljA.12|) for the monomials F{x) — x^ , x^, x^ , . . . we deduce, by induction on 771, that the 
coefficients Qm,\{0) in M'^^'^ are equal to the Taylor coefficients a„i of L^(5) for m = 2, 3, 4, . . .. □ 

Appendix B. From exponentials to quasi-polynomials: alternative approach 

One can also deduce the Euler Maclaurin formula for a polytope directly from an Euler Mclaurin formula 
for a cone. This approach is a bit longer than the approaches taken in Section |31 Here we outline this 
approach and include a lemma that may be of independent interest. 

The exact Euler Maclaurin formula for an exponential function on a non-singular convex polyhedral cone 
is this. Let F be a vector space with a lattice Vz, and let V£ c V* be the dual lattice. Let m,. . . ,Un be 
primitive elements of V£ which form a basis for V* , and let ai, . . . , a„ 6 be the dual basis. Take any 
Ai, . . . , A„ G Z and let w = — ^j'^j- Consider the finite abelian group 

r = V^^Vspan^lu,}. 

Consider the cone 

C„ = {a; I {uj,x) + \j >0 , j = 1, . . . , 71} = {w + ^ t^a^ \ tj > , j = 1, . . . , n}, 
and its expansions, given by 

Cy{h) = {x I {uj,x) + Xj + hj > , j — 1, . . . ,n} 

for h near 0. Let 



fix) 



where ^ G satisfies, for each j = 1, . . . ,n, 

(a) Re((^, aj)) < 0, and 

(b) I aj) I < 2tt\ {y, Uj) \ for aU y S V£ such that (y, aj) ^ 0. 

(The set of ^'s that satisfy (b) is a neighborhood of the origin in VJ.) Then 

n d 



(B.l) II = ^e2^«(7,") -Q 



h=0 



/• 
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This formula follows directly from I4.7|l . (|4.11|) . and H5.18|l . 

By applying this formula to the polarized cones Cj, that occur in the polar decomposition (Section (SJ, 
together with some bookkeeping, one deduces the exact Euler Maclaurin formula on a polytope, 

where ^ G V* is sufficiently small and is "polarizing", i.e., belongs to the complement of a finite union of 
(real!) hyperplanes through the origin. 

One would like to obtain a similar formula for polynomial functions by taking the derivatives of (|B.2|) 
with respect to ^ and taking the limit as ^ — > 0, (or, alternatively, by comparing the coefficients in the Taylor 
expansions in ^ of the left and right hand sides of (|B.2|) '). For this one needs to show that the infinite order 
differential operator TdA(gf^, • ■ • , gf^)i applied to /^(^^^ e^^'^^ dx, commutes with derivatives and limits with 
respect to ^. This follows from the following lemma, which may be of independent interest. 

Lemma B.l. Consider the exponential function 

(B.3) /(e,:,) = e<«--> 

where x and ^ e V^. Let bi, . . . ,bd be positive numbers, and let Tj\(Si, . . . , Sd) be a formal power series 
that converges on the multi-disk 

(B.4) \S,\<h , i = h...,d. 

Then the series 

(B.5) Ta(^,...,^)/ fi^,x)dx 

aril dtid JA(h) 

is absolutely convergent whenever satisfies the inequalities 

(B.6) I {£,, ai^v) I < bi for all v E Vcrt(A) and all i G 1^ , 

and this convergence is uniform on compact subsets of the domain (|B.6|) . 



Proof. By general properties of power series, Ta(5'i, . . . , Sd) is absolutely convergent on the multi-disk HB.4|I . 
and this convergence is uniform on any strictly smaller multi-disk. Let 

(B.7) TA{Si,...,Sd) = ^ ^ CpM Yl Si^^"* + terms that involve other monomials . 

nii non-negative integers 

The terms that involve other monomials make zero contribution to (|B.5|I . by H6.3|) . 

Because the series (jB.7|) is absolutely convergent on the multi-disk ljB.41) . so is the series 



E E n ^r- 

Ui non-ncgativc integers 

for each 5 G Z^^, where we set Cp^fii = if < for some i G Ip- For n = {ni)i(zjp and m = {mi)i^ip we 
write n < m to mean n.i < rui for all i E If, and, for such n, we write |m — n| = ^ TOj — ti^. Then the series 

/ \ 



EE 



|C'F,m| 

n such that 
n<m and 
\|m— n|<dim Y / 



is also absolutely convergent on the multi-disk ljB.4p . 

Let K be any compact subset of the set of ^'s that satisfy (|B.6|I . Choose positive numbers A and 6^ such 
that < 6'j < bi and < A < 1 and such that 

(B.8) I {£_, I < A6- for all v e Vert(A) and all i e ly 
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for every ^ E K. To prove the lemma, we will show that the series (|B.5p is dominated on i^T by a multiple 
of the converging positive series 

/ \ 



EE 

F n 



|C'_F,m| 

m such that 
n<m and 
\|m— n|<dim V / 



ieiF 



For each face of A and each i G Ip we choose an element ai^p G to be equal to ai^y for an arbitrary 
vertex w e F. Then the elements ai^p of V have the following properties. 

(1) 

(B.9) {ui,a.,p) = \ 

(2) By (|R8|) . for every ^eK, 

I cii.F) I < for each face F and index i G 7^?. 
The Stokes formula (|6.1|) and part (1) of (|B.9|) imply that for each face of A and each i € Ip 

(B.IO) 



/ / = - E '<^l^^^,P) / / Da,,ff- 

The exponential function IjB.Sfl satisfies -Da/ = a) / for any a G V. Combining these facts with 1().4|) 

, d 



and HB.10|I . we get 



(B.ll) 



_d_ 
dhi 



,OtiP) 



dhi 



A{h) 



I-{S,.a,.p)l f 

A{h) 



I /--E(-' 

Applying Hie/ aT (|B.11() . using the fact that the g^'s commute, and applying Ht).2|) . we get 

/ w . V, . . \ I t ...... ^ t 



dh^ \Tp\ Jpf^h) 



[ / = - E {uua^^p)-^ [ / - (C, ai,F) / , 

■Inh) ,«„.Kt,w |iF|7i.(,.) 



Z such that 
£;:=_Fncri 
satisfies IS^EC-F 



Iterating this formula we get, for any ii, . . . ,ik G Ip, 



li,...,ls such that l<ji < ...<js </c 

_B^:=_Fn(Tijn...n<T;^ 

satisfy 
F=_Eo2...2-Es7^0 

i2-l 



Let B > 1 be such that 



■ • ■ {ui^^ai-^^E, 

\{ui,ai^E)\ < B forain,i, and E'. 
By this and l|B.8|) . the term the right hand side of (|B.12|) that corresponds to some fixed h, . . . ,ls and some 
fixed ji, . . . , is is bounded by 



(B.13) 



j = l,...,k 

J0{ii,---,is} 
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where Eg ^ F Ci ai-^ H . . . Ci cri^ . Let Bi be strictly greater than B'^X '^jp^ |/e -^1 for aU < s < dim V and 
aU faces E of A. Then, for h near 0, the bound (|B.13p is less than or equal to 



(B.14) 



Si A* 



n 



j=i,...,k 



Let rrii be the number of times that i occurs in (ii, . . . ,ik)- Then the left hand side of ljB.12|) can be 
rewritten as 



nil- 



/• 



Let rii be the number of times that i occurs among ij for j G {1, 
can be rewritten as 



, fc}\{ji, . . . , js}. Then the bound HB.14|) 



(B.15) 



Denote rfi = (m^, i G Ip) and n — {rii, i e Ip). Then rii < rrii for all i, which we write as n < m, and 
^(rrii — Ui) < s < diniF, which we write as \rn — n\ < dimi^. Then we can further bound (jB.lSp by the 
following number which depends on m and not on n: 



(B.16) 



BiX max 

n such that 
n<m and 
|m— n|<dim F 



This bound was for the summand of ljB.12(l which corresponds to a fixed choice of li,...,ls and of 
ji, . . . , js- The number of possible choices of h,. . . ,ls is bounded by d**, which is bounded by d'^™^, and, 
further, by d*^™^. The number of possible choices for ji, . . . ,js is (J), which is bounded by fc*^™^, and, 
further, by fc*^""^. It follows that the right hand side of (|B.12|I is bounded by 



(B.17) 



^dim V ^f^irn V 



^lA max 

n such that 
n<rh and 
|m— n| <dim F 



Since 



and since the maximum among positive numbers is bounded by their sum, there 



exists i?2 such that HB.17|I is bounded by 
(B.18) B2 



^ E 

n such that iGIf 

n<.rh and 
|m— n|<dim F 



To conclude, we found B2 such that 



A(/0 



<B2 j: Uib^) 

n such that iGIf 

n<.rh and 
|m— n|<dim F 



The series (IB. 51) can be written as 



d 



l+nii 



/• 



p - J A(h) 

By what we have shown, for ^ in the compact set A', this series is dominated by the positive series 



F m 



n such that iGlp 

n<rh and 
|m~n| <dim F 
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which can also be re-written as 



( 



F n 



I c'f,??! I 



m such that 
ra>n and 
\|m — n| <dim V 



As we have shown, this series is convergent. □ 

References 

[AW] J. Agapito and J. Weitsman, The weighted Euler Maclaurin formula for a simple integral polytope, Asian J. Math. 9 
(2005), 199-212. 

[B] A. Barvinok, A Course in Convexity, Graduate Studies in Math. 54, Amer. Math Soc, 2002. 

[BP] A. Barvinok and J. E. Pommersheim, An algorithmic theory of lattice points in polyhedra, New perspectives in algebraic 
combinatorics (Berkeley, CA, 1996-97), 91-147, Math. Sci. Res. Inst. Publ., 38, Cambridge Univ. Press, Cambridge, 
1999. 

[BDR] M. Beck, R. Diaz, and S. Robins, The Frobenius problem, rational polytopes, and Fourier- Dedekind sums, J. Number 

Theory 96 (2002), no. 1, 1-21. 
[Be] J. BernoulU, Ars conjectandi Basel (1713). 

[Bo] N. Bourbaki, Functions d'Une Variable Reele, Chapitre VI: Developpments Tayloriens Generalises, Formule Summa- 
toire d'Euler-Maclaurin. 

[BVl] M. Brion and M. Vergne, Lattice points in simple polytopes. Jour. Amer. Math. Soc. 10 no. 2 (1997), 371—392. 
[BV2] M. Brion and M. Vergne, Residue formulae, vector partition functions and lattice points in rational polytopes, Jour. 

Amer. Math. Soc. 10 no. 4 (1997), 797-833. 
[CSl] S. E. Cappell and J. L. Shaneson, Cenera of algebraic varieties and counting lattice points. Bull. A. M. S. 30 (1994), 

62-69. 

[CS2] S. E. CappcU and J. L. Shaneson, Euler- Maclaurin expansions for lattices above dimension one, C. R. Acad. Sci. Paris 

Sr. I Math. 321 (1995), 885-890. 
[CS3] S. E. Cappell and J. L. Shaneson, unpublished. 

[dCP] C. de Concini and C. Procesi, The algebra of the box-spline, arXiv:math.NA/0602019l 

[DR] R. Diaz and S. Robins, The Ehrhart Polynomial of a Lattice Polytope, Ann. Math. (2) 145 (1997), no. 3, 503-518, 
and Erratum: "The Ehrhart polynomial of a lattice polytope" , Ann. Math. (2) 146 (1997), no. 1, 237. 

[Du] J. J. Duistermaat, The heat kernel Lefschetz fixed point formula for the spin-c Dirac operator. Progress in Nonlinear 
Differential Equations and their Applications 18, Birkhauser Boston, Inc., Boston, MA, 1996. 

[Eu] L. Euler, Commentarii Acad. Petrop. 6 (1732-3) and 8 (1736). 

[F] Paul Filliman, The volume of duals and sections of polytopes, Mathcmatika 39 (1992) no. 1, 67-80. 

[GV] A. N. Varchenko and I. M. Gelfand, Heaviside functions of a configuration of hyperplanes. (Russian) Funktsional. 

Anal, i Prilozhen. 21 (1987), no. 4, 1-18, 96. EngUsh translation: Functional Anal. Appl. 21 (1987), no. 4, 255-270. 
[GGK] V. L. Ginzburg, V. Guillemin, and Y. Karshon, Moment maps, cobordisms, and Hamiltonian group actions, to appear 

in the American Mathematical Society University Lecture Series. 
[Gr] B. Griinbaum, Convex Polytopes, Second Edition, Graduate Texts in Mathematics 221, Springer- Verlag, New York 

2003. 

[GLS] V. Guillemin, E. Lerman, and S. Sternberg, Symplectic Fibrations and Multiplicity diagrams Cambridge University 
Press (1996) 

[Gu] V. Guillemin, Riemann-Roch for toric orbifolds, J. Diff. Geom. 45 (1997), 53-73. 

[Ha] C. Haase, Polar decomposition and Brion's theorem, math. CO/0312095 Contemporary Mathematics 374 (2005), 
91-99. 

[H] G. H. Hardy, Divergent Series, Second Edition, 1991 (unaltered), Chelsea Pub. Col, New- York, N.Y. 

[HM] A. Hattori and M. Masuda, Theory of multi-fans, Osaka J. Math. 40 (2003), 1-68. 

[KSWl] Y. Karshon, S. Sternberg, and J. Weitsman, The Euler Maclaurin formula for simple integral polytopes, Proc. Nat. 

Acad. Sci. 100 no. 2 (2003), 426-433. 
[KSW2] Y. Karshon, S. Sternberg, and J. Weitsman, Euler Maclaurin with remainder for a simple integral polytope, 

|arXiv:math"C O /0307125 

[KKl] J. M. Kantor and A. G. Khovanskii, Integral points in convex polyhedra, combinatorial Riemann-Roch and generalized 

MacLaurin formulae, Inst. Hautes Etudes Sci. preprint. 
[KK2] J. M. Kantor and A. G. Khovanskii, Une application du theoreme de Riemann-Roch combinatoire au polyndme 

d'Ehrhart des polytopes entiers de R'^, C. R. Acad. SciParis Ser. I Math. 317 (1993), no. 5, 501-507. 
[KW] P. J. Kelly and M. L. Weiss, Geometry and Convexity: a study of mathematical methods. Pure and applied mathematics, 

John Wiley and Sons, Wiley-Interscience, New York, 1979. 
[Khl] A. G. Khovanski, Newton polyhedra and toroidal varieties, Funktsional. Anal, i Prilozhen. 11 (1977), no. 4, 56-64. 

EngUsh translation: Func. Anal. Appl. 11 (1977), no. 4, 289-296 (1978). 
[Kh2] A. G. Khovanski, Newton polyhedra and the genus of complete intersections, Funktsional. Anal, i Prilozhen. 12 (1978), 

no. 1, 51-61. English translation: Func. Anal. Appl. 12 (1978), no. 1, 38-46 (1978). 



EXACT EULER MACLAURIN FORMULAS FOR SIMPLE LATTICE POLYTOPES 



35 



[KPl] A. G. Khovanskii and A. V. Pukhlikov, Finitely additive measures of virtual polytopes, Algebra and Analysis 4 (1992), 

161-185; translation in St. Petersburg Math. J. 4 (1993), no. 2, 337-356. 
[KP2] A. G. Khovanskii and A. V. Pukhlikov, A Riemann-Roch theorem for integrals and sums of quasipolynomials over 

virtual polytopes, Algebra and Analysis 4 (1992), 188-216, translation in St. Petersburg Math. J. (1993), no. 4, 789-812. 
[Kii] K. Knopp Theory and application of infinite series, Translated from the Second German Edition by Miss R. C. Young 

L.es Sc., Blackic and Son Limited, London and Glasgow (1928), Chapter XIV. 
[L] J. Lawrence, Polytope volume computation. Math. Comp. 57 (1991), no. 195, 259-271. 

[Ma] C. Maclaurin, A treatise of Fluxions Edinburgh (1742). 

[Md] L. J. Mordell, Lattice points in a tetrahedron and generalized Dedekind sums, J. Indian Math. Soc. (N.S.) 15 (1951), 
41-46. 

[Mo] R. Morelli, Pick's theorem and the Todd clas of a toric variety. Adv. Math. 100 (1993), no. 2, 183-231. 
[Pi] G. Pick, Geometrisches zur Zahlentheorie, Sitzenber. Lotos (Prague) 19 (1899), 311-319. 
[Poi] S.D. Poisson, Memoires Acad, sciences Inst. Prance 6 (1823). 

[Po] J. E. Pommersheim, Toric varieties, lattice points and Dedekind sums, Math. Ann. 295 (1993), 1-24. 

[S] J. Shancson, Characteristic Classes, Lattice Points, and Euler-MacLaurin Formulae, Proceedings of the International 

Congress of Mathematicians, Zurich, 1994. Basel: Birkhauser Verlag, 1995. 
[SV] A. Szenes and M. Vergne, Residue formulae for vector partitions and Euler-MacLaurin sums. Adv. in Appl. Math. 30 

(2003), no. 1-2, 295-342. 

[Ti] E. C. I'itchmarsh, The Theory of Functions, Second Edition, Oxford University Press, 1950. 

[V] A. N. Varchenko, Combinatorics and topology of the arrangement of affine hyperplanes in the real space. (Russian) 
Punktsional. Anal, i Prilozhen. 21 (1987), no. 1, 11-22. English translation: Functional Anal. Appl. 21 (1987), no. 1, 
9-19. 

[Vel] M. Vergne, Residue formulae for Verlinde sums, and for the number of integral points in convex rational polytopes. 

European women in mathematics (Malta,2001), 225-285, World Sci. Publishing, River Edge, NJ 2003. 
[Ve2] M. Vergne, private communication. 

[Wi] W. Wirtinger, Einege Anwendugen Der Euler- Maclaurin'schen Summenformel, insbesondere auf eine Aufgabe von 

Abel, Acta Mathematica 26 (1902), 255 - 271. 
[Z] G. M. Ziegler, Lectures on Polytopes, Graduate Texts in Mathematics 152, Springer- Verlag, New York, 1995. 

Department of Mathematics, The University of Toronto, Toronto, Ontario M5S 3G3, Canada 
E-mail address: karshon@matli.toronto.edu 

Department of Mathematics, Harvard University, Cambridge, MA 02138, USA 
E-mail address: shlomoSmath.harvard.edu 



Department of Mathematics, University of California, Santa Cruz, CA 95064, USA 
E-mail address: weitsmanSmath. UCSC.EDU 



